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Abstract 



This paper considers the quantum collapse of infinitesimally thin dust 
shells in 2 + 1 gravity. In 2 + 1 gravity a shell is no longer a sphere, 
but a ring of matter. The classical equation of motion of such shells in 
terms of variables defined on the shell has been considered by Peleg and 
Steif [l] , using the 2 + 1 version of the original formulation of Israel [2] , 
and Crisostomo and Olea [3], using canonical methods. The minisuper- 
space quantum problem can be reduced to that of a harmonic oscillator 
in terms of the curvature radius of the shell, which allows us to use well- 
known methods to find the motion of coherent wave packets that give 
the quantum collapse of the shell. Classically, as the radius of the shell 
falls below a certain point, a horizon forms. In the quantum problem one 
can define various quantities that give "indications" of horizon formation. 
Without a proper definition of a "horizon" in quantum gravity, these can 
be nothing but indications. 



It is a great pleasure to be able to contribute to this special volume in honor of 
Octavio Obregon. One of us [M.R.] has known Octavio for more than 30 years 
and has collaborated with him on many fruitful projects. We both want to wish 
him many more productive years and continued success. 



1 Introduction 

This article describes a toy model of a toy model. There has been some interest 
over the years in minisuperspace quantization of thin shells as models of the full 
quantum collapse of more complicated objects. One of most the vexing problems 
in this scenario has been whether an event horizon will form and consequently 
some sort of "quantum black hole," or whether a shell of non-interacting parti- 
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cles will simply collapse to a point where the uncertainty principle will provide 
a "repulsive force" and the shell will reexpand into the same universe. 

Simple considerations can answer this question up to a point. It seems 
obvious that for large enough masses one might expect a reasonably peaked wave 
function centered on a radius near to but outside the classical Schwarzschild 
radius of the shell would, as it moves toward zero radius, maintain enough 
coherence to pass beyond the Schwarzschild radius almost in its entirety and 
a "black hole" would form with high probability. For small masses one might 
expect that rapid spreading would overwhelm the coherence of the wave function 
and the shell would reexpand into our own universe with probability essentially 
(or exactly) one. 

In full relativistic quantum gravity these simple ideas are fraught with diffi- 
culties. The most serious of these are: 

1) The quantum problem becomes unphysical for masses much above the 
Planck mass. While the minisuperspace approach has frozen out all radiative 
modes of the gravitational field and one is insisting on a single-particle interpre- 
tation of the shell problem, the wave function in the Schrodinger picture can still 
become pathological at a point where one might expect graviton production to 
begin. Haji'cek [1] has given a sufficient but not necessary condition that shows 
that we might expect problems above a couple of Planck masses. In [S], a qual- 
itative argument was presented that used Compton wavelength considerations 
to give a similar bound. 

2) There are very serious technical problems in the formulation of the prob- 
lem. Many authors have used the full ADM method to construct a Hamiltonian 
for the system, where they have chosen an internal time for the system in order 
to have a true Hamiltonian and a Schrodinger equation that allows the study of 
the time evolution of the quantum system. The choice of an internal time leads 
to the problem of quantum formulations that are not unitarily equivalent. An- 
other problem with these Hamiltonian formulations is that they often require an 
ad hoc choice of a Hamiltonian in terms of variables defined on the shell itself. 
Several of these Hamiltonians have been given by various authors [1], [5]. A 
Hamiltonian due to Hajfcek and Kuchaf [T] has the advantage of being defined 
by a coherent procedure with no ad hoc choices, and is formulated in terms of 
foliations of spacetime by timelike surfaces. I will discuss this Hamiltonian in 
more detail below. All of these Hamiltonians have limits on the mass of the 
shell of a few Planck masses. One formulation that does not seem to have a 
mass limit is based on the Wheeler-DeWitt equation corresponding to one of 
the Hamiltonians in [6j. 

3) Many of the Hamiltonians are quite complicated and there is no real 
chance of finding analytic solutions to the Schrodinger equations of these models. 
Numerical solutions have been presented by a group consisting of A. Corichi, 
G. Cruz, A. Minzoni. M. Rosenbaum and M. Ryan of the UNAM, N. Smyth of 
the University of Edinburgh, and T. Vukasinac of the University of Michoacan 

In order to study a simpler problem (another toy model of a toy model), 
Ryan considered the quantum collapse of a shell in Newtonian gravity. It 
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is possible to derive a classical equation of motion for the radius of the shell in 
terms of just the radius, R{T), and the shell mass, M. This problem has the 
advantages that there is only one time and that the equation is the same as that 
of a particle falling radially toward a point particle of mass M located at r = 0. 
It is easy to quantize this system, and the final Schrodinger equation is the 
analogue of the s-state hydrogen-atom equation. While exact analytic solutions 
exist, one has to use scattering states, and the integrals needed to form wave 
packets are not tabulated, so only approximate solutions were given in [H]. 

The idea of the present article is to study the quantum collapse problem in 
2 + 1 gravity, where one can address some of the difficulties mentioned above 
in a context where some of the problems mentioned above do not exist. The 
quantum problem is unambiguous for all masses, so there is no problem of wave 
function pathologies. As will be shown below, one possible Hamiltonian for this 
problem has the form of that of a harmonic oscillator. This will allow us to 
find simple analytic solutions that can be used to illustrate the development of 
the quantum collapse of the shell, and will allow us to investigate the problem 
of horizon formation in terms of analytic constructs. Our toy model of a toy 
model will help us to understand exactly what is happening in the numerical 
solutions that will be presented in [S]. In the final section of the article we will 
mention possible extensions of the 2 + 1 problem that might be the subject of 
future research. 

The plan of the rest of the article is as follows. Section 2 will be a brief resume 
of the literature on the general relativistic minisuperspace problem. Section 3 
will set up the classical 2 + 1 problem, Section 4 will consider the quantum prob- 
lem and present solutions, while Section 5 will be conclusions and suggestions 
for future study. 

2 The collapse of thin shells in relativistic quan- 
tum gravity 

The study of the quantum collapse of dust shells is about a decade old [TD]. 
The classical collapse problem is fairly straightforward, and can in principle be 
solved exactly. One assumes a 5-function massive (or null) shell where Birkhoff 's 
theorem tells us that outside the shell the metric is Schwarzschild and the metric 
inside the shell is Minkowski. The Israel junction conditions can be used to 
derive the equation for the evolution of the shell in terms of intrinsic variables 
on the shell itself, the proper time, r, of an observer riding on the shell and the 
curvature radius, R{t), of the shell that he would measure. The equation for 
the motion of the shell becomes [2] 



where M is the rest mass of the shell (a constant of motion) and M is the 
Schwarzschild mass of the exterior metric. It is straightforward to define x = 




(1) 
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MR, V = dx/dr and M/M as our Hamiltonian. Assuming V — V{P) and 
solving dH/dP = V{P) for the "momentum" P, we find that V = sinh~\P) 
and our Hamiltonian becomes 

H = coshP , (2) 

2a;' ^ ^ 

where m = A4/Mpi, Mpi the Planck mass, a Hamiltonian given by Hajfcek 
[3]. Unfortunately, this is not the only Hamiltonian that gives the equation 
of motion ([T]), and we are left with the problem of defining an "appropriate" 
Hamiltonian for the problem. A number of Hamiltonians for different choices 
time, including the time of an observer at the center of the shell where space is 
flat, and a Wheeler-DeWitt equation identical to that for a relativistic charged 
particle radially falling in a Coulomb potential are given by Haji'cek, Kay and 
Kuchaf [6]. 

Kuchaf and Hajfcek [7], dissatisfied with such ad hoc Hamiltonians, have 
managed to construct a Hamiltonian for collapsing dust shells that comes di- 
rectly from an ADM reduction of the Hilbert-plus-matter action. The problem 
with this approach is that, as spacetime quantities, the matter variables are pro- 
portional to (5[i? — Ro{t)], where Ro{t) is the position of the shell as a function 
of shell proper time. It is virtually impossible to reduce the time derivatives of 
these delta functions to reasonable variables in the matter Lagrangian that can 
give a shell Hamiltonian that describes the motion purely in terms of canonical 
variables on the shell. Kuchaf and Haji'cek used an ingenious method based on 
the fact that the ADM reduction by a 3 -I- 1 foliation is not restricted to foliation 
by spacelike surfaces, but works just as well for foliations by timelike surfaces. 
Using this approach and a formulation of the dust fluid velocity in terms of 
velocity potentials, they define a new Hamiltonian. The cost of this consistent 
formulation is a very complicated Hamiltonian, 

1/2 



H = -V2i?|^I-^-^I-^cosh|j R>2M, (3) 

1/2 



/ j\ T 2M P \ 

H=-V2R[i J Isinh- 0<R<2M. (4) 

\ R y R R I 

Almost all of the Hamiltonians that have been given (except for the Wheeler- 
DeWitt equation of [6]) seem to have mass limits beyond which the wave func- 
tions become pathological. These Hamiltonians are all so complicated that it 
seems impossible to find analytic solutions to assist in the interpretation. In 
[S], Corichi et al. have given a series of numerical solutions that give the evolu- 
tion of wave functions that are initially sharply peaked over a radius near the 
classical horizon, R = 2M, for the coshP Hamiltonian ([2]) and the Hamiltonian 
given by Eqs. ^S^. All of these solutions show evolution of the peak toward 
R = with a bounce caused by the boundary conditions at i? = with the 
appearance of interference fringes as well as a rapid spread of the wave packet. 
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The Kuchaf-Hajicek Hamiltonian has wave functions similar to those of but 
with many rapid oscillations superimposed. 

Since these Hamiltonians are self-adjoint, unitary evolution implies that a 
peak formed from scattering states will always rebound to i? = cx). One can 
ask whether this behavior means that all quantum collapse of this sort implies 
a rebound into our own universe. Since r is proper time on the shell and R is 
also a shell variable, such questions can only be answered by knowing the global 
quantum spacetime surrounding the shell. In any case, the scenario of the shell 
observer is that he sees (begging questions of quantum measurement and the 
reduction of the wave packet) the shell collapse to some point near R = 0, 
where uncertainty principle effects change the classical equations of motion and 
the shell rebounds (actually, a shell where the particles do not interact directly 
with one another "passes through itself" and reexpands, that is, each radially 
infalling particle passes through R ~ and the azimuthal angle 9 jumps from 
the initial 9o to 6o + tt). Even if the shell has collapsed below its classical 
horizon, in finite proper time it will again be above the horizon and traveling 
toward R = oo. This quasi-classical scenario is not surprising. In proper time, 
a classical shell that manages to avoid forming a curvature singularity at i? = 
would behave in this way, but as R becomes greater than 2M the shell would 
be in a universe past our temporal infinity («+), or in "another universe." 

This quasi-classical scenario is what one might expect to see for a large mass 
where the quantum fluctuations would be small compared to the mass and the 
evolution of the wave packet would be coherent long enough for the shell to 
collapse past its horizon and the shell would emerge from the horizon into a 
new universe. However, if the wave packet spreads sufficiently so that the width 
is greater than the classical horizon radius, we can see that a horizon might 
never form, and the shell would reexpand into our own universe. 

The problem of horizon formation in the quantum system is very difficult. 
Event horizons are global features and one has to try to define a global feature in 
a fluctuating manifold. Of course, this quantum manifold must be constructed 
in terms of the full minisuperspace canonical quantum gravity of the shell- 
metric system. In the shell case we tend to use some kind of approximation 
to construct the spacetime metric. Kuchaf [TT] argues that for the simple shell 
minisuperspace we may just replace the the shell mass (Schwarzschild mass) and 
the shell radius in the metric outside the shell by the corresponding operators 
to make a "metric operator." The problem with this metric operator is that it 
is a function of shell proper time, and studies of the metric close to the shell [12] 
cannot tell us whether a true event horizon (tied to observer time at infinity) 
forms. 

Other approximations are under study [Sj. The simplest calculation is to 
calculate < R{t) > and the uncertainty AR = \J < (R— < R >)2 > and check 

whether AR becomes very large as < i? > becomes small so that < R > does not 
fall below the classical horizon and AR is larger than the classical horizon, which 
can be taken as an indication of the non-formation of a horizon. In [S] , numerical 
evaluations of these two quantities will be presented for the Hamiltonians ^ 
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and ([3J1]), and they suggest that no horizon forms for small masses. Another 
possibility (to be considered in [9]) would be to take M\tjj{R, t)P to be a classical 
density p{R,t) and calculate the classical metric due to a classical fluid with 
this mass distribution and see whether a horizon forms. Note that M should be 
either the rest mass or the Schwarzschild mass. It is not yet clear which. There 
are technical problems with this calculation. We have to calculate the metric 
from a density that is given in terms of a solution of the Schrodinger equation 
for our Hamiltonian, and there is no guarantee that this density can be made to 
obey the equation T'^''^ — for our fluid. Another possibility considered in [3] 
is that of the "metric operator" mentioned above, which was extended to the 
whole manifold outside the shell and used to define a "quantum" stress-energy 
tensor. 

The rotationless 2 + 1 problem has some advantages over the 3 + 1 problem. 
The Hamiltonian can be constructed fairly easily, and, as will be shown, has the 
form of a harmonic oscillator. The Schrodinger equation for this Hamiltonian 
has well-known analytic solutions. The expectation value of R and its uncer- 
tainty can, in principle, be calculated analytically. The analogue of the other 
calculation using p{R, r) is much simpler than in the 3 + 1 case. The classical 
horizon is easily found. In the next section these ideas will be considered. 



3 The classical problem of shells in 2 + 1 gravity 

The first element we need for this problem is an equation for the radius of the 
shell. This problem has been studied in detail by Peleg and Steif [T], using the 
2 + 1 version of the original formulation of Israel [5] , and Crisostomo and Olea 
[3], using canonical methods. We will review the calculation of Peleg and Steif 
(using the parametrization of the metric of Crisostomo and Olea) . 

The calculation we will use follows that of Israel who, as mentioned above, 
studied the collapse of a shell in 3 + 1 gravity, represented by a delta-function 
sphere of dust of radius R{t). As we have mentioned, in 2 + 1 gravity the shell 
is a circle, i.e. a ring of matter, also of radius R{t). The metric of spacetime 
will be written in circular coordinates, where flat space is represented by the 
metric 

ds^ ^ -dt^ + dr^ + r^d9^. (5) 
The equation for a circle R{t) is 

(3)r = i?(r), (3)t = i(r). (6) 

We will use the notation i,j = 1, 2, 3 and A,B = 1, 2. We will now need a set 
of coordinates on the circle, which we will take to be S^a — {ttO), and we 
will define two three- vectors on the circle by e^^^ = dx'^/d^A, where are the 
coordinates t, r, 9 given above. We have 

eU = {^,R,0), (7) 
e(e) = (0,0,1). (8) 



6 



The metric on the circle is (c?s^ — gABd£,^dS,^) 

ds^ ^ -dT^ + R^{T)d9^. (9) 

We now have to calculate the three-dimensional metric due to the ring of 
matter. We have to solve the 2 + 1 Einstein equations (necessarily with cosmo- 
logical constant A to avoid completely flat solutions), 

R^J~lg^JR=-A9^J, (10) 

inside and outside the circle. Since there is a Birkhoff theorem in 2 + 1 gravity, 
these metrics in vacuum will be static or stationary. We will study the static 
case, where the matter has no angular momentum. In this case the metric has 
the form 

ds^ = -f{r)df + JJ^dr^ + r^de^, (H) 
and the well-known solution is 

/ = S-Ar^ (12) 

B a constant. Since A has units of inverse length, we will write, as is common, 
A = ±l/£^. The final form of the static, circularly-symmetric metric is 

ds^ = - f 1 - 2M T dt^ + j-dr^ + r^dO^, (13) 

where we have taken B, following Criostomo and Olea, to be 1 — 2M, M a 
Schwarzschild mass. Notice that this is not the choice of Bafiados, Teitelboim 
and Zanelli (BTZ) [13j, who take B = — A/btz, but we choose 1 — 2M following 
Cristo somo and Olea because we will want the metric inside the ring to be the 2 
-I- 1 AdS metric. These two choices of B correspond to a shift in the zero of the 
mass parameter [M], [15]. The relation between these two masses (and others) 
and the ADM mass is discussed in the Appendix. This choice requires that 
the mass term be negative in order to have a 2 -I- 1 black hole with a horizon. 
Notice also that we must have AI > 1/2 in order to have a horizon. As in the 3 
-|- 1 case, we expect that outside the ring we will have a black hole metric with 
some "Schwarzschild mass" M, that is, 

ds' - - f 1 - 2Af + ^) dt^ + ^dr^ + r^de^, (14) 

V (l-2M + f^) 

and inside the ring we will have a metric that is as near as possible to flat space 
(but with a cosmological constant). This will be p4)) with M — 0, that is the 2 
+ 1 AdS metric, 

ds^ = -(l + ^] dt^ + —^-^dr^ + r^de\ (15) 
V (1 + iJ) 
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With these preHminaries we can follow the steps of Israel's derivation of the 
equation of motion of a shell in 3 + 1 gravity, a manner similar to that used 
by Peleg and Steif [Tj. For any three-dimensional vector Aa, we can define 
components of on the surface of the circle as 

AA = A,el^y A'=A^elAy (16) 

We will also need the vector n' , the unit normal to the circle R{t), with n^rii = 
+ 1. We can now use the standard equation 



dA 



■7^ = A^'.AelB) - A^KBAn\ (17) 



where Kab is the extrinsic curvature of the ring, to arrive at the equation 
or 

n,u\^jU^ = -u^u^KABf, (19) 

where ± means the quantity calculated for r > R{t) (+) and r < R{t) (- 
). In order to calculate the left-hand-side of we will need expressions 
for and rii. Following Israel, = R, = Q and UiU^ = — 1 gives us 

= X = \J f + R^ / f ■ There are many ways to find n^, and we obtain nt = — ^, 
rir = X and ne = 0. We will also need the definition, 

lAB = K+s - K^B, (20) 

and the Lanczos relation, 

lAB - gABl = SttS'ab, (21) 

where qab are the components of the induced metric on the surface in terms of 
r and 0, j — g^^jAB and Sab is the surface stress-energy tensor. In our case, 
we will be interested in a dust shell, so we will take Sab = ctuaub, cf the rest 
mass density on the ring. 

Following Crisostomo and Olea, the Einstein equation for GijU^n^ becomes 

^R + K^bK±^ -Ki = -2G,jn'n^ = -Ag,jn'n^ = -A. (22) 

Subtracting the minus equation from the plus equation, and defining Kab = 
h^^AB + ^As)' arrive at 

jabK"^"" - kl - 0, (23) 

where K ~ qabK^^ . The Lanczos relation gives us SabK^^ — or, since 
Sab = a-u"^u^ , 

u^u^Kab = 0. (24) 
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Israel uses Eg. p^ . calculating the left-hand-side directly, and, doing the 
same for our case, we find 



{niu\jU^)^ = 



I 



(25) 



Now, summing the plus and minus versions, and using u^u^ Kab = 0, we find 

1 1 



0. 



(26) 



Using ^AB = ^T^iSAB ~ QabS) and calculating ^abu^u^ — u^u^ {Sab 
qabS = (t — (t = 0, we find that 



R 



R 



1 



1 



= 0, 



(27) 



which is consistent. Unfortunately, Israel uses the 3-1-1 equivalent of (P?]) . 
where the right-hand-side is not zero, to find a second-order equation of motion 
for R(t) in terms of the rest mass density, but our equations no longer contain 
any reference to this variable. However, we can use, following Peleg and Steif [T], 
the Lanczos relation, Eq. (l2T|) . directly. If we impose the condition R+R/P = 0, 
we have that [riiU^.^u')^ are both zero, so from Eq. (fT9|) . u^u^ K^g = 0. From 
we find that u"^ = (liO)j so K^^ — 0. Now, from the Lanczos 

and grr ~ ^U so 

(28) 



UA 



relation we have that Jtt — grrl = SttStt, but •frT 
Stt = f, and finally, 

7 = Sttct. 

Since 7 = g^^jAB = -Itt + (l/i?^)7ee = {l/R^)^ee, we have 

76/6/ = STiR^a. 



(29) 



\r=R' 



We can now calculate Kge directly from its definition, Kgg = rig^g — Tgg 
Calculating the necessary r*j, and using the fact that ng^g — 0, we find that 
Kgg = -fXR, SO 

Jgg^{fX)-R-ifX)+R. (30) 

From S^.^ = 0, we can show that a = m/27ri?, m the total rest mass of the 
ring. 

If we now consider the equation of motion of _R, R = —R/P, we see that it 
has a first integral, 

(31) 



E = 

2 



2^2- 



Using (fx) R — {fX)^R = 4mR, and inserting pTjl in this relation, we find 
VTT2E - Vl - 2M + 2E = 4m, (32) 
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which can be solved for E as 



E^— —^ V2m^--^ — [ h8m - -. (33) 

2 4to)2 2 2 32 V"i y 2 ^ ' 



We would like to use Equation (|3T|) to construct a Hamiltonian formulation 
of the problem. This equation is simply the energy equation for a harmonic 
oscillator, so it is obvious that if we take E as our Hamiltonian we can define a 
momentum Pr as R, and we have 

(34) 

Hamilton's equations for this Hamiltonian are equivalent to the equation of 
motion R^-R/f-. 



4 Quantum collapse in 2 + 1 gravity 

4.1 Schrodinger equation and evolution of wave packets 

Once we have a Hamiltonian in terms of the ring variables R and r, we can try 
to construct a Schrodinger equation for the problem. Before we do this, we have 
to study the units of our variables. If we write the classical relation 

we see that the quantities on the right-hand-side have no units, so the left-hand- 
side must also be dimensionless, so M and m must be dimensionless. We would 
like to multiply the true mass by some combination of G and c that will give 
a dimensionless mass. We have to be careful about the law of gravitation in 
two space dimensions. In general. Gauss's law would give us, in d dimensions, 
F = GMm/r'^^^ , and in two space dimensions F = GMm/r. This means that 
G would have units of [M]-'^[LY[T]-^ , and G/c^ has units of 1/[M], and GM/c^ 
is dimensionless. If we consider the metric function /(r) = 1— 2M-|-i?2/£2^ which 
has to be dimensionless, this means that M in conventional units is GM/c^. It 
is usual (see Ref. [TS]) to define a "Planck mass" Mpi = c^/G and a "Planck 
length" Lpi = TiG/c^ , even though there is no particular reason to identify this 
mass and length with the 3 + 1 Planck mass and Planck length, especially if we 
have no knowledge of what G should be in two space dimensions. 

In our case, M means M/Mpi and m means m/M-p\. In conventional units, 
our equation for R becomes 

Our "Hamiltonian" should have the units of energy for a conventional quantum 
Hamiltonian, so we should multiply Eq. (|35p above by Mpic^, our i?- momentum 
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Pr becomes Mp\R , and with 
Mpic2 ^ ^ ' 



H 



and 



4m / 



p2 



Mpi 



Mpic2 



^'^Pl 2 d2 



(37) 



(38) 



2Mpi 2 
with uo = c/£. 

We want to quantize this system with H replaced by an operator H. If we 
look at ((37)) . we see that the right-hand-side must become an operator, so the 
left-hand-side must also be an operator. One way to achieve this is to take AI to 
be a q-number, and m and Mpi c-numbers. Of course, there is no special reason 
to make this choice, but we do so to make contact with previous work (sec, for 
example, [H]). If we consider Eq. ([2]), this choice is motivated by the fact that 
the Schrodinger equation for ^ is similar to the hydrogen atom Schrodinger 
equation, with m playing the role of e^, and it is usual in the hydrogen atom to 
take as a c-number rather than a q-number. 

We can now take the wave function of the system to be 5* = ^'ji/4'(i?, r), 
with ^m{M) an approximate eigenstate of M with eigenvalue Mg. An ex- 
act eigenstate of M of this type would be 6{M — Mq), but to avoid problems 
with the integral of the square of a delta function we will assume that '^m 
is an extremely sharply peaked wave function centered on M ~ Mq. In this 
case, our Schrodinger equation becomes (M^^'a/ k, M^^m-, and realizing Pr 
as 



-ifid/dR) 



Mo 
4m 



2mc^ 



/ m 



\Mpi 



Mpic? 



ih 



5*(i?, r) 



cj^i?^*(i?,r). 



(39) 



2Afpi ai?2 

One curious fact about this equation is that, if we define = exp{—iET /h)ip{R) 
then the energy eigenvalues are 



En = (n + -)hujQ = (n + -) — , 



(40) 



and since E is given by the first line of ([39]) , there is a discrete relation between 
the Schwarzschild mass, Afo, and the rest mass, m, 

2 



Mq 2m 

— - H 

4™ Afpi 



-l = (2n + l)g = (2n + l)^|^. (41) 



We will now return to units where G = c = ?i = 1, Afg now meaning Mq/Mpi, 
m now meaning m/Mpi, and £ meaning l/Lp\. If we solve for Mq in terms of 
m and £, we find 



Mq = 8m^ 



1 



1 



Am? 2™^ 



1 



1 



n'^+2 



(42) 
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We now want to write our Schrodinger equation in terms of dimensionless 
variables. For a moment we will return to conventional units, and define the 
following dimensionless variables. We define y = R/^/lLpi and T — ct/1. The 
Schrodinger equation now becomes 



* dT dy^ 
which has eigensolutions 

VI/ = p-»(«+i)T 1 



y2^(7r)i/4 



y'^iy,T) (43) 



'""'^H.M, (44) 



Hn Hermite polynomials. 

Now, in order to make contact with previous work, we would like to study 
the evolution of a wave packet sharply peaked around a value of y ~ yo, some 
distance outside the point where a classical horizon would form if the radius 
of the shell were to fall below Rh = £y/2MQ — 1 and follow its movement as 
the packet falls toward R = 0. In previous work it was necessary to solve this 
problem numerically, since in the 3+1 case analytic solutions to the Schrodinger 
equation did not exist [5] [9], and in the Newtonian case, [8], while analytic 
solutions existed, it was impossible to sum the expansions for the relevant wave 
functions to give an analytic expression for the collapsing wave packet. In the 
2 + 1 case, however, we can give an exact analytic expression for the wave 
packet as a coherent harmonic-oscillator state. Even though the eigensolutions 
are nothing but harmonic oscillator wave functions, the radial variable R cannot 
be negative. In order to keep this from happening we will take the potential 
to be that of a half oscillator with an infinitely hard wall at i? = 0. This 
potential is shown in Figure 1. This will mean that ^(0,t) will always be 
zero. This boundary condition can be enforced by only expanding in odd-n 
harmonic oscillator eigenfunctions. We will want to begin with a difference of 
two Gaussian states, one peaked around y = +yo and the other around y =^ —yo, 
so that their sum at y = is zero. This state is 

0) = a[e-i'^y-y°^^ - e-*(2/+ao)=^]^ (45) 

only valid for y > 0, with a a normalization constant. Since this is a sum of 
two Gaussian states, we can use standard techniques to construct the difference 
between two coherent Gaussian states with (|45l) as an initial condition. The 
result is 

*(y,T) = ae-'-°^/2e*?^o-"2^°^/4[exp(-i{(y- yoCOSCc>oT)2})exp(-*2/2/osinc^or) - 

- exp(-i{(y + yo cosujot)"^}) exp(iyyo smujor)], (46) 
which is zero at y = for all r. The normalization a is easily found to be 
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We would like to connect the variables that describe the wave function with 
the radius of the classical horizon, Rh, shown in Fig. 1. We will make the 
following definitions, yo = XVl\/2Mo - 1, y = wVi\/2MQ - 1, T = cr/i. If we 
take \/£ — 1 1 y/2Mo — 1 , the unnormalized probability density, p = 5'*^'/q;^, 

p ^ ^-{w-XcosTf _|_ ^{w+XcosTf _ 

-2e-"''e-^' ^ cos(2A?i; sin T), (47) 

is shown, for A = 3, as a function of w for various values of T in Figures 2-7. We 
begin with a Gaussian packet at T = that collapses toward w = 0, developing 
interference fringes as it reaches a minimum for some w < X, then rebounds 
toward w = X again. This pattern then repeats forever. In order to compare 
our results with earlier work, we will only be interested in one cycle of this 
pattern. 

4.2 The formation of a horizon 

In previous articles the quantization of the shell collapse was used to study the 
possibility of the formation of a horizon in the quantum collapse. As mentioned 
in Sec. 2, this concept has many difficulties. An event horizon is a global 
construct and it has no local definition. This means that in quantum gravity 
one would have to return to the starting point and try to define what a "quantum 
horizon" might be. Once this definition has been decided upon, one must try to 
find out if some collapse process will result in the formation of such a horizon, 
with the result being a probability of horizon formation. Theories of quantum 
gravity in their present state are far from being able to give us this result, 
so in shell collapse some articles [5], [9], [8] have tried to give an estimate of 
horizon formation by finding out if a sharply peaked wave packet, during its 
collapse toward R ~ 0, will fall, in some sense, below the classical horizon 
radius, Rh- In some sense, because the packet will usually spread and basically 
will never lie entirely below R = Rh- One could use the integral of \I'*vI' from 
R = to R — Rh, which is a number less than one, as the probability of 
horizon formation. Another possibility would be to use the operator M in the 
expression Rh = — 1 to define a "horizon operator" that could be used 

to define the probability of horizon formation. In previous work it was decided 
to use < i? > as a function of r, which will fall from i?o to a minimum and then 
begin to increase again. If this minimum is below the classical horizon, one can 
say that a horizon forms and if not, not. This is a yes- no answer instead of a 
probability, but it is a quick estimate. We will use this concept below. 

It is not difficult to calculate < R > (t) from the wave packet given in 
but the result is a complicated function that contains the error function and 
Dawson's function, so trying to find the minimum of < i? > (r) by finding the 
point where d < R > / dr — requires the solution of a transcendental algebraic 
equation, making it difficult to give an analytic expression for the point where 
a horizon would form. Instead we will use the fact that, as can be seen from 
Fig. 5, that at T = 7r/2, the point where the packet begins to rebound, the 
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peak nearest i? = is high and narrow. We can use the position of this peak 
as a parameter to tell us whether a horizon forms or not. If the position of the 
peak is below Rh, a horizon forms, and if not, not. In previous work it was 
found that for large masses the shell collapse was so rapid that the wave packet 
fell below Rh so quickly that quantum mechanics did not allow it to rebound 
before that point, while for small masses the rebound occurred for R > Rh- In 
order to calculate the position of the maximum of p nearest to w = 0, we write 
p as 

piy,T) ^ 2e"^""^"™"'^[cosh(2yyocosr) - cos(2|/yo sinT)]. (48) 
The packet is closest to y = when T = tt/2, which gives 



p{y, Tr/2) = 2e-y [1 - cos(2yyo)] = 46"^' sin^iyyo)- 



(49) 



This function has zeros at y = mr/yQ, and either peaks or valleys (or inflection 
points) at dp{y,'K/2)/dy ~ 0, or 



8e y sin(yyo)[-ysin(?/?/o) +2/oCOs(yyo)] = 0. 



(50) 



An obvious solution to this equation is y = nir/yQ for the zeros. Another 
possibilty is y — yc, where 

yc ^ yocot{ycyo), (51) 

which are peaks (it is easy to show that (Pp/dy^ < 0). If we define yc = 
jVeV2Mo - 1, we have 

7 = Acot[7A^(2Mo - 1)]. 



for A = 3 we can solve for Mq as 



1 



(52) 



(53) 



We now have two relations, ((53)) and (Mo/4m+2m)^ — 1 = E. The minimum 
for Mq [the right-hand-side of decreases monotonically with 7] is when 
7 = 3, since we want yc l£ yo- The numerical value of this minimum depends on 
£, and there is no simple reason to choose any special value for £. We can use 
the relation between Mq, m and E to find another value for Mq as a function 
of i, and equating these two, we get a numerical value for i. There are various 
ways to define Mo as a function of £ from the energy relation. If we consider 
the classical equation 



Mo 
Am 



2m 



r 

e ' 



(54) 



and take j/o to be where Py = 0, then (Mo/4m 
One solution of this quadratic equation is 



2m)2 - 1 = y2£ = isMo - 9. 



Mo = 8m^ 



17- 



288 J 



(55) 
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Mo 



This is real for 18m^ > 1/18, which gives a minimum for Mq equal to Mq 
17/18, and for 7 = 3, 

1 / j_ ^\ 17 

or £ = 7r/32 « 0.098, and Mo(7 = 3) = 17/18 « 0.944. 

If we use the expectation value of H as our energy, we have < H > 
{Mo/Am + 2mY ~ 1, and calculate < H > to he 



18' 



(56) 



<H>- 



1 



(57) 



in principle we can solve (using i/q = 9£[2Mq — 1]) for Mg, but the transcendental 
equation is difficult to solve for Mq as a function of m and i. If we were able to 
assume that yo were sufficiently large that the the factor (1 + e~^o)(l — e~^o)~^ 
were one, we could solve the resulting quadratic for Mq easily. However, for 
the minimum Mq from ([53]) with 7 = 3, we find that Uq = 7r/4, and our factor 
is almost equal to three. We can find an approximate solution by inserting 
Vo = 7r/4 in the factor and write the quadratic equation 



Mn 



4m 



1 



-^ + 2to -l = - + (18Mo-9)(l 



,-7r/4N-l 



(58) 



which gives, using /3 = (1 + e-'"/'*)(l - er"/'^)-'^ = coth(7r/8) 
8m^ 



2.68 



Mo 



18/3- 




(18/?)2 - 36/3 - 



1 

Am? 



(9/3-1) 



which has real solutions for 



^^,^(9/3-1)^1/. 



(18/3)2 - 36/3 ' 



(59) 



(60) 



which has a minimum at = [(9/3 — 1) — l/£]/[(18/3)2 — 36/3], and inserting 
this expression in (j59p . and equating the result to (j53p with 7 = 3, we find 
^(l/9/3-l)[7r/3/4+(18/3-l)/(9/3-l)] « 0.180. This gives us Mo (7 = 3) = 0.742. 

For the classical approximation we now have Mo (7 = 1) = 2.62, while for 
< H > we have Mo (7 = 1) = 1.66. Since < H > should give an answer closer 
to the correct one, we can take it as giving the best value of £. This means that 
for 

0.742 < Mo < 1.66, (61) 

no horizon forms, while for Mq > 1.66 one does. For example, for 7 — 1/2, 
Mo = 2.80. 

From the Appendix, we can give these masses in terms of A/btz , and for 



1.48 < Mbtz < 3.32, 
no horizon forms, while for Mbtz > 3.32 one does. 



(62) 
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Of course, as has been mentioned above, the rebound of the wave packet, 
once it has passed the horizon does not mean that the shell returns through the 
horizon into the same spacetime where it began. In spite of the fact that in 
terms of proper time the shell exits the horizon in a finite time, in terms of the 
time, t, of an observer at infinity this exit occurs at time after t = oo, or into 
"another universe." 

5 Conclusions 

The 2 + 1 calculation given above has the advantage of being the only model 
so far that has given analytic solutions to the quantum collapse problem. This 
is a great advantage in understanding the origin of features in the numerical 
solutions given before in Refs. [S], [5], [H]j especially the the existence of the 
interference fringes that were noticed in Refs. [S] and [S]. The results for horizon 
formation are at least consistent with those of previous work. 

Another possible model problem would be to study a simple relativistic 
model of gravity such as the Nordstrom theory. In a such a scalar theory the 
reduction to a Hamiltonian formulation is relatively straightforward, but one of 
us (L. O.) (see Ref. [16]) has shown that the Hamiltonian is more complicated 
than that of ordinary 3 + 1 gravity. He has also shown that linearized gravity 
gives the same Hamiltonian as full gravity. 

Still another possibility would be to carry out the program suggested in Ref. 
[9] in the 2 + 1 case, that is, to try to define an average metric by writing the 
metric outside the ring in terms of coordinates ^ (where r = R{t) + ^) and 
t = T. In this metric we can replace R{t) by i?, to form a metric operator gij 
and calculate < gij > r), where, for example, ggg = {R + ^)^, so 

<gee>=<R^>{T) + 2<R{T)>^ + e, (63) 

and use the classical Einstein tensor, G^, calculated from < gij > as SirTij 
to define a classical stress-energy tensor Tij{^,T) (which automatically satisfies 
T^''.j = 0). From this Tij, we can define, for example, a shell density as a 
function of ^ and t, which could be plotted. 

Appendix 

The mass parameter M that appears in the black hole metric ITil) can be as- 
sociated with the ADM mass of the black hole. The ADM mass for the metric 
with B ^ ~M is defined as [14]. [15] 

Madm fdrde V^g^'^R « ^^A/btz, (64) 

Idtt J r 8 

where ^'^^gAB is the two-metric on surfaces of constant t, and is the Ricci 
scalar on these surfaces. The factor £/r is related to proper time at r ^ oo [15] 
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and is usually ignored, so we find that 



Mbtz = x;— Madm- (65) 

We can choose our mass unit to be any multiple of Afpi. In and [151 they 
take it to be 8Mpi while in [IJ it is taken to be simply Mpi. 
If we calculate A/adm for our metric, we find that 

MADM = T^-(M-i) + A (66) 

where A is a constant taken to be (Mpi/4)r/£ in [M] and [15] that represents 
the zero of the ADM energy, Madmc^. In older work [13], BTZ chose A = 
[Mpxj A)r 1 1 + (Mpi/8)^/r as we have in this article. The difference lies in the 
metric for M = 0. In our case, M = is AdS, while for later BTZ work M — 
represents 

ds^ ^ --prdt^ + -^dr^ + r'^de'^. (67) 

In our case, the choice of reference mass equal to Afpi leads to a relation 
between our M and Mbtz, 

M = hiBTZ. (68) 
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Figure 1: The harmonic oscillator potential for our problem with M-p\ijj'^l2 = 1. 
The dashed line shows a typical position of the classical horizon radius, Rh- 
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Figure 2: The probability density for the wave packet given in Eq. (|47|) for 
r = (the initial state). 
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Figure 5: The probability density for the wave packet given in Eq. ()47|) for 
T= 1.57, «7r/2. 
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Figure 6: The probability density for the wave packet given in Eq. (|47|) for 
T = 2.08. 
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Figure 7: The probability density for the wave packet given in Eq. (|47|) for 
T = 2.6. 
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Figure 8: The probability density for the wave packet given in Eq. (|47|) for 
T = 3.14, WTT. 
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